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It is popular to probe the structure of the QCD vacuum indirectly by studying individual fermion eigenmodes, 
because this provides a natural way to filter out UV fluctuations. The double-peaking in the distribution of the 
local chiral orientation parameter (X) has been offered as evidence, by some, in support of a particular model of the 
vacuum. Here we caution that the V-distribution peaking varies significantly with various versions of the definition 
of X. Furthermore, each distribution varies little from that resulting from a random reshuffling of the left-handed 
(and independently the right-handed) fields, which destroys any QCD-induced left-right correlation; that is, the 
double-peaking is mostly a phase-space effect. We propose a new universal definition of the X parameter whose 
distribution is uniform for randomly reshuffled fields. Any deviations from uniformity for actual data can then 
be directly attributable to QCD-induced dynamics. We find that the familiar double peak disappears. 



1. Introduction 

The use of eigenmodes to probe the structure 
of the QCD vacuum is fairly old, but the re- 
cent resurgence of interest is the result of the fo- 
cus on local properties P and the use of chiral 
fcrmions |2I3| . The space-time chiral structure 
of the eigenmodes reflects the space-time topo- 
logical structure of the underlying gauge fields. 
Our work 0] has involved a detailed study of the 
space-time geometric structure and we have ar- 
gued that this is inconsistent with the instanton 
liquid model (ILM). Many others |S] have studied 
the behavior of the low lying eigenmodes and have 
argued in favor of the ILM, but their conclusion 
relies solely on the interpretation of the observed 
double-peaked behavior of the probability distri- 
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bution of the local chiral orientation parameter. 
To help resolve this, here we revisit the local chi- 
ral orientation parameter. 

2. Local Chirality Parameter 

The local chiral parameter, X, measures the 
tendency of a field at a space-time point to 
be left or right handed. It was originally 
defined in terms of the polar coordinate 9 = 
tan -1 l^l/IV^I- Many definitions of X are pos- 
sible (including several from the literature [S]): 



A« = -tan 
X^ = Itan 



1^ 



1 = X 



-i M 



L|2 



R\2 



X^ = -tan 

L|2 



X 



BBC 



|V£|Va 

IVfl 2 



l#l 2 + hA 



B\2 



1 



2 



A = +1(— 1) for a purely left (right) handed field. 
However, as we see in Fig. 1, the double peaking 
depends strongly on the definition of A. We need 
a universal measure which will produce a constant 
distribution for random, uncorrelated data. 

V=16 4 a=0.082 fm f=000-002 pair 01-02 




Figure 1. Alternate definitions of X. 

3. A Test: Randomly Reshuffling Fields 

To determine how much of the double peaking 
in the A-distribution is due to dynamics of the 
QCD vacuum, we make the following test: Com- 
pute the X distribution from ip L (x) and tp R (x). 
Then randomly reshuffle the fields: ipra,a( x ) = 
ip L (x lan ), where x lan is a random permutation. 
Use an independent random reshuffle for ip R . 
Now any QCD-dynamically-induced correlation 
between ip^kn and ipf an at site x is destroyed. Re- 
compute the A-distribution. We see in Fig. 2 that 
the randomized A distribution looks very similar. 
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Figure 2. Original A-distribution and the A- 
distribution after random reshuffling. 



To understand how there can be so little 
difference between original and randomized A- 
distributions, consider the marginal probability 
p(ip L ) — J dip R p(ip L , ip R ) °f the actual (non- 
randomized) Monte Carlo data (Fig. 3). The key 
is that the marginal distribution has a long tail. 
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Figure 3. Marginal probability distribution for 
\ip L \ or \4> R \. 

To construct an A distribution, one first cuts 
the data, selecting the most intense points, i.e. 
the top 2% (or 5% or 10%) of the density tpty = 
ip L ^ip L + %jj R ^ip R . For -0 L and ij) R chosen ran- 
domly from this distribution, it is likely that ei- 
ther ip L or ip R is large (to make the cut). Given 
this pre-selected data, and one of these large com- 
ponents (say ip L ) in the tail of the distribution, 
it is very likely that the randomly chosen ifj R is 
much less than the selected ip L (since the bulk of 
the marginal distribution is in the peak at smaller 
values). But such a pair (large tfj L and smaller 
ip R , or vice versa) result in a value of A much 
different than 0, leading to the "double peak". 

4. A New Universal Definition of A 

The "A-distribution" is a histogram of the 
probability p(A) of a given value of A. With 
generic coordinates qi = \ip R \ and <?2 = |^ ,L |, 



dq 1 dq 2 p(qi, q 2 )S (X - X(q 1 , q 2 )) 



But if the components q\ and q 2 were uncorre- 
lated (supposing there were no dynamics), then 
one would obtain a different histogram, 



dq\dq 2 p(qi)p(q 2 )S (A - X(q 1 ,q 2 )) 
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Clearly, to expose the true correlation we need to 
"subtract" the phase space background. Consider 
the cumulative distribution function P(X) = 
J^dX' p(X'). Then P G [0,1]. For uncorre- 
lated components, what is the probability density 
of obtaining a particular value of P? 

p(P) = ^ dX p(X)S (p- J ^dy p(y)j 
= 1 

That is, this probability density is uniform for 
uncorrelated data, which is what we seek! We 
can rescale and rename this quantity X{X) = 
2P{X) - 1 so that X G [-1, +1] as for X. X(X) 
is the improved quantity against which we should 
plot our distributions. It will be constant for un- 
correlated data. Deviations from uniformity will 
appear for correlated data. 



5. Results 
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Figure 4. For various definitions of X: Origi- 
nal A"-distribution. X-distribution after random 
reshuffle. "Difference" (plot to. new variable X). 

In Fig. 4 for each of four definitions of X, we 
plot the original data and its randomization ver- 
sus X, and their "difference" (in thick black) ver- 
sus the new variable X. Compare to the baseline 
of no correlation (thin black line) at 0.02 (50 his- 
togram bins). Note that the new improved def- 
inition X is universal; that is, the black curves 
are exactly the same. Thus the "Double-Peak" is 
exposed as a phase-space effect, which is removed 
with the new X. 



6. Summary 

Fermion eigenmodes (which filter UV fluctua- 
tions) are useful probes of the structure of the 
QCD vacuum. Observation of a double peak- 
ing in the local chirality parameter {"X distri- 
bution" ) is a necessary but not sufficient condi- 
tion for the viability of the instanton liquid model 
(ILM) 4 . The double-peaking depends strongly 
on the definition of X. Moreover, after the fields 
have been randomly reshuffled (thus destroying 
any QCD-dynamics correlations), the X distri- 
bution looks very similar to the original distribu- 
tion! That is, most of the double peaking is due 
to phase space and does not reflect the dynamics 
of QCD. Accordingly, we have constructed a new 
universal measure ( U X" ) of local chirality against 
which one should plot probability distributions. 
It is uniform for uncorrelated data. Deviations 
from uniformity appear for correlated data. The 
"double-peak" , exposed as mostly a phase-space 
effect, is flattened with the new improved X, leav- 
ing little evidence in support of the ILM. 
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